Quantum fidelity and quantum phase transitions in matrix product states 
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Matrix product states, a key ingredient of numerical algorithms widely employed in the simulation of quantum 
spin chains, provide an intriguing tool for quantum phase transition engineering. At critical values of the control 
parameters on which their constituent matrices depend, singularities in the expectation values of certain observ- 
ables can appear, in spite of the analyticity of the ground state energy. For this class of generalized quantum 
phase transitions we test the validity of the recently introduced fidelity approach, where the overlap modulus 
of ground states corresponding to slightly different parameters is considered. We discuss several examples, 
successfully identifying all the present transitions. We also study the finite size scaling of fidelity derivatives, 
pointing out its relevance in extracting critical exponents. 
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I. INTRODUCTION 



The study of many-body quantum systems is at the same 
time a fascinating and challenging field of research. This is 
due to the richness of inherently complex phenomena aris- 
ing in the presence of a large number of interacting particles, 
among which quantum phase transitions (QPTs) occupy a dis- 
tinguished position. These transitions take place at zero tem- 
perature and are driven by an external parameter, as for exam- 
ple the magnetic field in superconductors. 

Recently it has been shown that the quantum fidelity - the 
overlap modulus - of two finite size ground states correspond- 
ing to neighboring control parameters is a good indicator of 
QPTs. Indeed, the fidelity typically drops abruptly at criti- 
cal points, as a consequence of the dramatic state transforma- 
tion involved in a transition. This approach has been tested in 
various contexts where the ground state can be calculated ex- 
actly, either analytically or numerically. These include Dicke 
model, one-dimensional XY model in a transverse field, and 
general quadratic fermionic Hamiltonians'-^. 

In the present article we further pursue this approach by 
analyzing QPTs described by matrix product states (MPSs)^. 
These states are at the basis of efficient numerical methods 
used in the analysis of spin chain systems, as the density ma- 
trix renormalization group (DMRG) algorithm. When con- 
sidered dependent on a control parameter g they can give rise 
to generalized QPTs, i.e., transitions where some observable 
quantity presents a non-analytic behavior in spite of the regu- 
larity of the ground state energy"*. We show that the quantum 
fidelity of two neighboring (in terms of g) MPSs is an effec- 
tive tool not only in detecting the critical point gc but also in 
giving the correct scaling at gc- The success of the fidelity 
approach in analyzing MPS-QPTs further proves thegener- 
ality of the procedure, which, as stressed in Refs. qIj, does 
not require any a priori understanding of the structure (order 
parameter, correlation functions, topology, etc.) of the con- 
sidered system. It is also worth pointing out that this method 
seems to have some advantage with respect to other quantum 
information based techniques. For example, quantum phase 
transitions in MPSs do not give rise to the logarithmic diver- 
gence of the entropy of block entanglement observed in other 



systems , thereby ruling out this quantity as a reliable tran- 
sition indicator On the other hand, other entanglement mea- 
sures can be used, e.g., single site entanglement and concur- 
rence, whose derivatives often provide interesting information 
about criticality*'. We discuss the results obtained by these 
methods in one of the simple examples we analyze below. As 
it will be shown, also these latter measures turn to be less ef- 
fective than the fidelity approach. 

In the following we will first present a derivation of the 
overlap formula for general MPSs (Sec. discussing some 
general features of these systems, and then study in detail the 
MPS examples introduced in Ref. 4 (Sec.|III|. As in Refs. [IB 
the analysis will be carried out also in terms of fidelity deriva- 
tives, which are the most suitable tool to observe finite size 
scaling properties. In particular, we will explicitly show how 
to extract the critical exponent of the correlation length from 
these quantities (Subsec. IIII Al) . thus demonstrating the inde- 
pendence and the completeness of our approach. Finally, we 
will report on the entanglement analysis in Subsec. UlIBI 



II. FIDELITY FOR MATRIX PRODUCT STATES 



Suppose we have a closed spin chain with sites. Let d 
be the dimension of the Hilbert space at each site. The matrix 
product states are defined asi 



|.g) := 1^(5)) = ^7 



TT{A,,...A,^)\ii 
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Zl,...,!^— 

(1) 

where the Aj's, with j = 0, . . . , d — 1, are D x D matri- 
ces, D is the dimension of the bonds in the so-called valence 
bond picture, and Et.'.,^«=o |Tr (A,, . . . )|2 a 

normalization factor. Since our goal is to explore quantum 
phase transitions in such states, we assume that the matrices 
Aj = Aj{g) depend on one or more parameters, generically 
denoted by g. The overlap of two MPSs corresponding to dif- 
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ferent g's is given by 

{9i\92) = W{g,)N{92)V'^ (2) 
X ^ Tr [Al (51) . . . (5i)]Tr [A,, (g,) . . . A,„ (52)] , 

where i = (ii. . . . , iw)- In order to simplify the formulae, we 
will work with the unnormalized overlap 



F{gi,g2) ■■= Vm9iWi92) {9i\92 



(3) 



and we obviously have M{ga) — Fi9a, 9a), a = 1,2. Using 
the identity Tr (A)Tr (B) = Tr (A (g) B), the overlap can be 
computed exactly for MPSs^ and we obtain 



F{g„g2)= Tr [i?^ (51 , 32 )] = ^ 1 ' ^2) 

k=l 

where are the eigenvalues of the x matrix 

1=0 



(4) 



(5) 



The latter is the generalization of the transfer operator 
Fi{g) := E{g,g), which we assume to be diagonalizable 
with eigenvalues Xk{g) ■— Xk{9,g)- Also, since A® B and 
B ^ A Sire isospectral, one has F{gi,g2) = F{g2,gi)*- 

The main object of our study will be the fidelity of two 
neighboring states 



T{g-S) :=|(.9-% + <5>| 



\F{g-5,g + S)\ 
^M{g - S)M{g + 5) 



(6) 



for a small variation 6 in the parameter space spanned by g. 
Defined this way, the fidelity depends not only on the param- 
eter g driving the QPT, but also on its variation 6. Expanding 
the fidelity in S we obtain J^{g; S) ~ I + djT{g: 5)\^^^5'^ /2 
(since J-{g\ 5) reaches its maximum for 5 — Q the first deriva- 
tive vanishes). Thus the rate of change of the fidelity close to a 
critical point is given by the second derivative of the fidelity''^ 



S{g) ■.^djT{g-5)\,^^ 



(7) 



and this is the relevant quantity we will study to determine the 
scaling at the critical point g = gc- Note that S{g) is con- 
nected to the ground state variation \\p'{g)) = dg\'4>{g)). For 
example, for real states one has the simple second order ex- 
pansion J-{g] 5) ~\ — 2(5^1 1 ^''(5) I p. The general expression 
of the second derivative has the formic 



S{g)^-^dA l^F{g,,g2\ 
Indeed, from T{g]Q) = 1 and dsT{g]5)\ 



(5=0 



(8) 



one has 



d}\nT{g-5)\,^^ = 5)1,^0 = S{g) and then 



F{g-6,g + 6)Fig + 5,g-5) 



F{g-S,g-S)F{g + S,g + S) 



<5=0 



{d's^d'g)lnF{g-S,g + S)\^^^ 



(9) 



from which Eq. ((H]) follows immediately. 

The general behavior of the fidelity for matrix product 
states in the thermodynamic limit (TDL) can be inferred from 
Eq. (IDl. Without loss of generality we can assume the eigen- 
values Afc of the generalized transfer operator E to be or- 
dered |Ai| > IA2I > ... > |Ai32|. ThenF(5i,52) = 

Ef=i (51,52) = Af [1 + Ef=2(Afc/Ai)^] and, provided 
|Ai| > |Afe| for k — 2, . . . , D^, in the thermodynamic limit 
N ^ 00 one finds F{gi,g2) ^ A;j^. Since the same holds 
for the normalization factors Af{g) = F{g, g), whenever the 
maximum eigenvalues Ai are non-degenerate, the large N be- 
havior of the fidelity must be of the form J-{g; 6) ~ a^, with 
< a < 1 by construction. Typically, a ^ 1 and the fidelity 
decays exponentially in the TDL^. If E\{g) exhibits a level 
crossing in the largest eigenvalue for a critical coupling gc, 
the degeneracy |Ai(5c)| = |A2(5c)| gives rise to a discontinu- 
ity of some expectation values in the TDL (called a general- 
ized QPT in Ref. 4). In this case the previous discussion has to 
be modified in order to include all the degenerate eigenvalues. 
In general, the vanishing of the fidelity is typically strongly 
enhanced at critical points. 

A similar analysis can be done for the second derivative 
S{g) of the fidelity. When a single eigenvalue Ai dominates, 
in the TDL one can replace F with Af in Eq. (HJ, thereby 
recovering the expected^ linear scaling cx N: 



S{g) ^ -m dg.dg, lnAi(5i,52)L 



(10) 



In contrast, at the critical point two or more eigenvalues of 
E{gi,g2) are equal in modulus. For example, at 5 = gc one 
can have |Ai(5c)| = |A2(5c)| > |Afe(5c)| for A; > 2. Then in 
the TDL 



Sigc) 



lim 

lim 

a^aa 



A9g,9g.,ln(Af + A 



V2^) 



-AN 



2 



4AfAf 
(Af + Ar 

Af 



dg^dg^lnXk 



,(11) 



where A1.2 — Ai. 2(51, 52)- Hence, at the transition one typ- 
ically recovers the divergence S{gc) oc N'^ already observed 
in Ref. 2. We also note that, by defining 9 :— ln(Ai/A2), 
the term in the second line of Eq. ( fTTT l can be rewritten in the 
compact form [N\dg-^9\/ cosh(A^6'/2)]^. This follows from 
the equality dg^ \n{Xi/X2)dg^ ln(Ai/A2) \dg, ln(Ai/A2)p, 
due to the relation Ai^2(.92,5i) — AJ 2(51,52)- Then, it is 
clear that this term vanishes for N ^ 00 unless '3i6 — 0, i.e., 
I All = IA2I. In the latter case, instead, assuming dg-^9 to be 
finite, the scaling cx N'^ is evident. Finally, we point out that 
the limit 5 5c in Eq. ( fTTT l is necessary to deal with possible 
singularities in the derivatives of the A^'s at the transition (see 
the examples below). 



III. EXAMPLES 

In the following we are going to study the explicit examples 
of MPSs described in Ref. 4. Without loss of generality we 
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FIG. 1: Fidelity J-{g; S) = \{g — 5\g + 5) \ as a. function of g for 
Example 1 with 5 = 10"^ and TV = 1000, 2000, 3000. Inset: close 
to the critical point, the fidelity oscillates for states belonging to dif- 
ferent phases (S = 10"^ N = lO"*). 



assume gi = g — S < g2 — g + S. 



A. Fidelity and second order MPS-QPTs 

Example 1. Take D = 2,d = 3 and {Aj} = 
{—Z,(T^,ga+}; in the following we will use the notation 
X, Z for the Pauli matrices, a± being the corresponding rais- 
ing and lowering operators. This one-parameter family of 
MPSs contains the ground state of spin-1 AKLT model for 
g = ±2 and has a critical point at gc = (see Ref. @|). The 
model is also significant because it has non-local string order. 
The parent Hamiltonian with two-body interaction is given by 



+2(4-g'){S!r + {g + 2f{S^St+,r + 
+g(.9 + 2){5f5f+i,5,-S,+i}], 



(12) 



where Si is the spin operator acting on the i-th site. The eigen- 
values of the transfer operator E = Z (S) Z + a- ^ a- + 
gig2<T+ (g) (7+ are (— 1, — 1, 1 ± ^/gigT) and Eq. (|4|l yields 



52) = (1 + VamY' + (1 - Vfflff^) +2(-i)'' (13) 

and M{g) = F{g, g) = {1 + g)^ + (1 - 5)^ + 2(-l)^. 

Let us analyze the behavior of the fidelity T{g; S) defined 
in Eq. Since !F{g; S) is evidently symmetric in both g and 
S, we can safely assume g > Q, S > Q. Three cases can then 
be distinguished, according to the behavior of the eigenvalues 
of the transfer operators: 

(i) g > S, where J^{g; S) ~ for ^ 1, with a = 
-S){l+g + 6) < 1; 
a/(l + 2(5)^/2for7V > 1, with 
a = 2 for TV even and a = 25^N{N - 1) for TV odd (this 
is the case where the largest eigenvalue of E{g — (5, 5 + (5) is 



{i + ^grzrpf/^^TTV- 

(ii) g = S, where ^{6; S) ^ 



FIG. 2: Finite size scaling for S{g). The second derivative S{g) 
is plotted for TV = 100, 200, . . . , 500 (logarithm is base 10). The 
peak at g = scales with TV accordingly to Eq. J16t . The dashed 
line corresponds to the asymptotic behavior given by Eq. l ll5t . 
The inset shows the data collapsing for the same curves (TV = 
100, 200, . . . , 500) when plotted in rescaled units. 



degenerate and the decaying behavior of the fidelity is slightly 
modified); 

(iii) g < S, where two of the eigenvalues of 
E{gi,g2) are complex and, for TV 1, the fidelity ex- 
hibits an oscillatory behavior driven by |cos(TV(p)|, (p — 
arctan y^P — g^, with a decaying envelope given by 2a^ , 

a = ^il + S^-g^)/[il + S 



itly, at the critical point g = gc 



g){l + 6 + g)] < 1. Explic- 
= the overlap becomes 



c|„ (l + ^^)^/^cos(TVy.o) + (-l)^ 



[{1 + S)^ + {l-S)^]/2+{-l)^ 

where (^0 — arctan (5. Then !F{0;6) = \{—S\S)\ (note that the 
overlap can be negative at this point). 

The fidelity behavior in the neighborhood of the critical 
point is shown in Fig. [T] for S — 10^"^ and various values 
of TV. It is evident that, although for a fixed S the fidelity van- 
ishes in the TDL for any g, the rate of decrease is much faster 
at the transition. This fidelity drop provides a useful finite 
size precursor of the quantum phase transition taking place in 
the infinite system. In the figure inset we zoom on the oscil- 
latory behavior in the interval g G {—6,5), where states on 
different sides of the transition are considered. Note that os- 
cillations start to be visible only for Nip > tt, i.e., TV > TVqsc, 
^osc = tt/ arctan \/ — g^ , which for .g = and 5^1 
reduces to TVqsc ~ tt/J. 

The second derivative S{g) of the fidelity can be calculated 
analytically using formula ([8]l. We find 



5(ff-0) = 



-TV(TV- 1) 



(TV even) 

-(TV-2)(A^-3)/3 (TV odd) 



(15) 



,(16) 
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where Eq. ( fTsT i corresponds to Eq. ( fTOb and gives the asymp- 
totic behavior for large N, while Eq. ( fT6] l is exact. As in Ref. 2 
the critical point is identified by a divergence cx of the sec- 
ond derivative S{g) of the fidelity. 

In Fig.|2]we plot log |S'((7)| for different values of N. In the 
inset of the figure the curves are plotted in rescaled units, giv- 
ing rise to data collapsing (lines for different values of N are 
practically indistinguishable in the inset, merging into a single 
thick line). This shows that S{g) /N"^ is basically a function of 
gN only. This is reminiscent of the usual scaling behavior of 
diverging observables in second order quantum phase transi- 
tions. Indeed, if an observable P of a one-dimensional system 
diverges algebraically in the TDL for some critical coupling 
gc, i.e.. Poo ~ Coo|g — .9c|^'', the finite size scaling Ansatz in 
the critical region reads Pn = N^^'^ Q{N\g—gcY'), where i^is 
the critical exponent of the bulk correlation length and Q is 
some function^. This is due to the fact that at criticality, once 
properly rescaled by some power of the size of the system, the 
observable is expected to depend only on the dimensionless 
ratio L/Cc between the system size L oc N and the corre- 
lation length. Since the latter diverges at the transition like 
l/\g — oils immediately finds -L/^c oc N\g — gd'^. 
In our case Pn = S/N = NQ{N\g\) with gc = 0, so that 
we obtain p = h' — 1. The result p = I agrees with Eq. ( fTsT i, 
which, for g ^ 0, yields Poo ^ — \g\~^- The value of h' 
can instead be checked by explicitly calculating the correla- 
tion length'^, where one finds 1^ = 1, confirming the finite 
size scaling estimate. 

The next examples of this subsection correspond to three- 
body interaction Hamiltonians, again taken from Ref. 4. Qual- 
itatively, they all feature the same behavior as Example 1: 
in the asymptotic limit \S{g)\/N diverges proportionally to 
l/\g\ at the critical point gc = 0, while the peak height of 
S{g = 0) scales with A^^, 

Example 2. Consider!) — d — 2andylo = (1— Z)/2+o-_, 
Aq = {1 + Z)/2 + .gcr+. The three-body Hamiltonian is W.2 
symmetric and has a critical point at gc = 0, where the state 
is a Greenberger-Horne-Zeilinger (GHZ) state, and reads 



H = Y, 2{g^-l)Z,Z,+i-{l+gfXM9-^fZ,X,+iZ, 



i+2 ■ 



(17) 

The generalized transfer operator E{gi,g2) has eigenvalues 
(0, 0, 1 ± y/gig2), very similar to Example 1. The function 
F{g — 6, g + 6) is then still symmetric in g and S. 

Since the largest eigenvalue Ai is the same as in Exam- 
ple 1, the TDL behavior of J-{g, S) and S{g) is unchanged. 
However, the parity dependent term (—1)^ is now absent. 
Then, the asymptotic behavior of S{g ^ 0) is still given 
by Eq. (flSl l. while the parity dependence of Eq. ( fT6l l is lost 
and one simply has S{g = 0) -2N{N - 1). The lat- 
ter results follows from Eq. (fTTT i. which applies exactly to 
this case even for finite A^. Indeed here only two eigenval- 
ues of E{gi,g2) are different from zero and one has the level 
crossing Ai(f; = 0) = \2{g = 0). Note that, due to the di- 
vergence of the double derivative dg^dg^ In Ai_2 calculated in 
di = 92 = 0> also the term in the last line of Eq. (fTTT i gives 
rise to a contribution cx A^^, in spite of its apparent linearity 
in A^ for |Ai| — IA2I. This explains the necessity of the limit 



g — > gc in Eq. ( fTTT i. 

Example 3. Consider Aq ^ X, Ai ^ ^{1 - Z)/2. The 
eigenvalues of E{gi,g2) are (±1, (^51^2 ± V5i52 + 4)/2). 
As before, the scaling of the system can be seen from the 
behavior of S{g) in the thermodynamic limit, driven by 
Ai (31,52) = (V^^i52 + \/5i52 +4)/2. We find 



S{g^O) ^ 
Sig = 0) = 



AN 



|5|(52 + 4)3/2 ' 
-A^V4 (A^ even) 



(18) 
(19) 



where Eq. ( fTSl l gives the large A^ behavior The parity de- 
pendence of Eq. iT% is caused by the negative eigenvalues of 
E{gi,g2), giving rise to terms oscillating like (—1)^. 

Example 4. Take Aq = a+, Ai = (t_ + y/g(]l + Z)/2. This 
MPS is the ground state of the following Hamiltonian: 

H = — g{Xi + XjZjj^i + ZjXj^i + ZiXi+iZi+2) 



+ (1 + 2g^)Zi — 2ZiZi+i — ZiZi+iZi+2 ■ 



(20) 



The eigenvalues of E{gi,g2) are (0, 0, (^5152 ± 
y/9192 + 4)/2). Hence, the asymptotic behavior of S{g) is 
again given by Eq. ( fTSl ), while S{0) ~ ~N'^/2 for A^ even 
and S{0) = -{N^ - l)/6 for A^ odd. Parity dependence 
arises from {^/gTg2 ~ V.9i.92 + 4)/2 = -1 for gi = .92 = 0. 



B. Fidelity, concurrence, and single site entanglement 

While all the above examples fit in the same picture, which 
is typical of the second order quantum phase transitions al- 
ready studied in Refs. 1.2, the next example, albeit trivial as 
matrix product state, features a different behavior, which will 
serve as a basis for some additional comments on the fidelity 
approach to QPTs. 

Example 5. Take D = d = 2, Aq = ^ Ai ^ 



9 











1+g 



The matrix E{gi,g2) has the eigenvalues 1 + gi, 



1 + .92, (1 + .91) (1 + .92), and 1 + 5"52- discussed in 
Sec. ini QPTs take place when two or more eigenvalues of 
^1(9) = E{9i 9) share the largest modulus. We then have to 
compare 1 + 5, (1 + 5)^, and 1 + g^". Clearly, for g < one 
has Ai [g) — l+g^"", while for <? > the dominant eigenvalue 
is determined by the ratio r{g) := (1 + 5^") /(I + g)^- The 
equation r{g) = 1 can be rewritten as g{g^"^^ — g ~ 2) — 0. 
Apart from the trivial solution gc = 0, for n > 2 a second 
solution g'c > 1 exists, so that in the region < g < g'c one 
has Xi{g) = (1 + g)^, while for g > g'c one obtains again 



Ai(g) = l 



For n = 1 the point g'^ can be considered at 



infinity and the latter level crossing never occurs. 

In order to evaluate the fidelity we have to substitute the 
above eigenvalues in Eq. (|6]l. For simplicity, here we restrict 
our analysis to the fidelity of ground states belonging to the 
same phase, avoiding to discuss the effects shown in the in- 
set of Fig. [T]for Example 1. In the phase corresponding to 
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FIG. 3: Fidelity for the model given in Example 5, with 5 = 10 ^ 
and iV = 100; n = 1 (top) and n = 2 (bottom). 



the interval g £ {0,g[,), where Xi{g) = 
Ai(5i,52) = (1 + g - S){1 + g + 5). 
delity in the TDL then results J-(g; S) ' 



(1 + g)^, one has 
The normalized fi- 

1, so that one has 



g + S)/^il+g- 5)^(1 + g + 6)T 
a constant phase in this region, due to the factorization of the 
largest eigenvalue of E{g—6, g+S). Outside this region such a 
factorization is absent and the fidelity vanishes exponentially 
with the size of the system. 

The asymptotic behavior of the second derivative S{g) can 
be calculated accordingly. In the interval g £ (0, g'^), where 
the eigenvalue (1 + g)^ dominates, one finds S{g) = in the 
TDL, while for g ^ [0,5^]' where the leading eigenvalue is 
given by 1 + g-^", Eq. ( fTOl i yields the formula 



Sig) 



-ANn^g^''-^/{\ 



(21) 



The latter function, which is evidently symmetric in g, has a 
minimum at g^m = -[{n - l)/(n + l)]!/^" (since |.gnii„| < 
g'^, the minimum position for g > is instead given by g^ it- 
self, as S{g) increases monotonically for g > g'^). Note that, 
while S{g) is continuous at g = gc = 0, where S'(O) — 0, 
at the second critical point g = g^ one has the discontinuity 

2)/[5c(l + 5c)^]' where the latter results is obtained by substi- 
tuting the defining relation r{g'^) = 1 into Eq. (I2TI 1. Further- 
more, exactly at the critical point g^ one has the superexten- 
sive scaling S{g'^) ~ -N'^[2n - 1 + [n - 1)3^7(1 + g'^)^. 

The finite size behavior of the fidelity is shown in Fig. [3] 
for N — 100 and 6 = 10""^. Due to the approximated relation 
T{g;6) ~ H-S'(5)(5^/2 one can there recognize also the scal- 
ing of S{g). In particular, in the lower panel where the case 
n = 2 is plotted, the minimum at gmin and the superextensive 
scaling at g'^{n = 2) ~ 1.521 are evident. 

For this simple example it is also possible to find a com- 
pact analytic expression for the ground state \g). Due 
to the commutativity [Aq, Ai] = 0, one has ^/J7\g) 



■ In, 



with k = X^jLi^i- Then 



by using Tr(A^) = 
for k 0, one finds 

VAf\9) = [ 



1 + (1 + .g)^ and Tt{A^-''A'i) = g 



(1 



(l+5)^]|0)«^ 

gr\or^ + {i- 



N 

fe=i 



\D 



N I 



(22) 



where \D 



N I 



J ' EiG/, \H---iN) with h ^ {i e 

Z^|X;^i«i = fc} is a Dicke state'" and \(f){g)) = (|0) + 

.g"|l))/-\/l + 3^" is a normalized single site state. Since the 
normalization is M = TrE^, depending on the dominating 
eigenvalue of Ei three cases are possible in the TDL: 

(i) Ai(.g) = (1 + g)-' > \2ig) and \g) ^ |0)«^; 

(ii) Ai(5) = 1 + > A2(ff) and \g) ^ 

(iii) Xi{g) = X2(g) and, in the non-trivial case g = g^, one 
hasi^lg) ^ (|0)®^ + |0)®^)/x/2. 

The above analysis can be summarized in the following ta- 
ble, where we recall the results for the largest transfer operator 
eigenvalue Ai ( g) and for the large N behavior of the function 
S{g) and the ground state | g): 





5e [0,.g^) 


g = g'c 


gi[o,g'c] 


Ai(.9) 


a + g? 


1 + 5'" = (1 + 5)' 


1 + 


S{g) 







oc N 


\g) 




m^'' + mg'c))^'')/V2 


10(5))^^ 



It is worth noting that the apparent drop of the fidelity at 
.gmin in the g < phase for n > 2 does not correspond 
to a QPT. Indeed, the scaling of S{gynin) does not present 
any peculiar behavior with respect to the rest of the g < 
phase. On the contrary, the scaling of S changes its nature 
at g = and g ~ g'^, highlighting the transitions, although 
in a very different way from the typical behavior observed in 
the previous second order QPTs. The transition at g = cor- 
responds to a discontinuity in dg {X) (Ref. 4), while that at 
g — g^ to a discontinuity of {X) itselfii^. Due to the permuta- 
tion symmetry arising from the commutativity of Aq, Ai, the 
correlation length cannot here be defined (the usual formula 

= 1 / In I Ai / A2 1 does not hold). In practice, the transitions 
of Example 5 appear to be of first order. 

For this example we also investigate concurrence and sin- 
gle site entanglement, mainly focusing on the critical point 
gc — 0. Of course, the fact that the state is factorized in the 
TDL trivially implies the vanishing of these quantities in this 
limit. We are however interested in the finite size scaling of 
the derivatives of these measures, in the spirit of Ref. 6. We 
are motivated by various aspects. First, as mentioned before, 
in MPS -QPTs the block entanglement entropy does not ex- 
hibit the logarithmic divergence observed in other quantum 
phase transitions, posing the question whether other entangle- 
ment properties usually signaling quantum criticality behave 
differently in this case. Second, we are interested in com- 
paring the effectiveness of the fidelity approach against other 
quantum information-theoretical methods. The simple nature 
of this example allows for a fully analytical treatment. 

Concurrence. We compute the concurrence for the re- 
duced density matrix of 2 qubits in the chain (the choice of 
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the qubits is unimportant, due to the permutation symmetry 
mentioned above). The density matrix of two spins p'^^ is 
obtained by tracing the initial p = \g){g\ over all the other 
spins in the chain. By exploiting the symmetries arising from 
the commutation relations [A^jAj] = and recalling that 

P''^'' Si,j,fc,i=o pfj',ki (fc^l is a 4 X 4 real symmetric ma- 



trix, one finds pf^,^i = p^^[^i = p^/jk- Together with the 
normalization Tr p*^^) — 1 this implies that only 5 entries of 
are independent, e.g., p^^l^^^, p^^l^^, p^^l^^, p^^l^^, p'^^^, 
significantly simplifying the calculation in the standard basis. 
Otherwise, one can use Eq. ( |22] |. easily obtaining 



.(2) 



-,(2) _ 



[(1 + g)^^ |00)(00| + (1 + 5^")^ |#) (23) 

+ (1 + .9)^(1 + ff'")(|OO)(00| + mm)]/^{9) ■ 



Having p'^'^\ the concurrence for two qubits is defined as"^ 
C = max{0, ^//ir - ^/pj - ,/JIi - ^/Jh), where p^s are 
the eigenvalues, in decreasing order, of the operator p^^^ {uy ® 
® <^v)- We finally have: 



C{9) 



2g^"|(i + gn 



(24) 



In the thermodynamic limit both C{g) and its first derivative 
vanish, so that the kind of transition signatures found, e.g., 
in Ref. ^ are absent here. However, for n = 1, the fourth 
derivative dgC\g^Q has a divergence in the TDL; for n > 2 
the singularity is present in even higher derivatives. 

Single site entanglement. This quantity is the von Neumann 
entropy of a single spin S = — Tr(p(^) logj p^^^), where p*^^^ 
is the reduced density matrix of a single site. Again from 
Eq. ( l22b one finds 



p'^^ = [(1+5) 



27V 



+(i+.9)"v/rT?^(io)(0i 



(25) 

l0>(O|)]/AA(g), 



whose eigenvalues are A± = (1 ± Vl-4detp(i))/2, with 
A+ = 1 — A_ and 



detp(i) 



g^"(i + gr^[(i + g ^") 

AA2(g) 



(26) 



The single site entanglement entropy S = ~\+ logj A+ — 
A_ log2 A_ vanishes in the thermodynamic limit at g = 0, 
together with its first derivative. For n ~ 1, one finds a di- 
vergence in 9g5|^^Q, while the divergence is shifted to higher 
derivatives for n > 2, similarly to the concurrence. The con- 
sidered divergence is however due to the functional form of 
the von Neumann entropy. Indeed, A_ ^ for g 0, giving 
rise to a singularity in the logarithm. 



We conclude this analysis by briefly discussing the TDL 
of these entanglement measures at the critical point g = 
g'^, where the state is not factorized (see the table above). 
The single site reduced density matrix in the TDL is simply 
P^'\9'c) = (|O)(O| + |</.)(0|)/2anddctp(i)(.gi) = [1-1/(1 + 
g',f]/A. Then \±{g[) = [1 ± 1/(1 + .g^)]/2 and S{g',) + 0. 
Note however that C{g'^) = in the TDL, reminiscent of the 
behavior of the GHZ state. This can be seen from the simple 
TDLofEq. (|23]l, namely, p(2) (5^) = (|OO)(OO| + |0(/))(00|)/2. 
IV. CONCLUSIONS 

In conclusion we have investigated the quantum fidelity of 
slightly different matrix product states dependent on a param- 
eter g in the context of quantum phase transitions. For generic 
MPSs, the overlap can be related analytically to the eigenval- 
ues of a generalized transfer operator, among which the largest 
in modulus plays a crucial role. If the latter is non-degenerate 
(in modulus), the fidelity typically exhibits an exponential de- 
cay in the thermodynamic limit and its second derivative is 
proportional to the size A'^ of the system. If instead more 
eigenvalues share the largest modulus, a quantum phase tran- 
sition can take place and the fidelity second derivative S{g) 
generally scales with A^^. We have demonstrated this behav- 
ior in the exhaustive analysis of some simple examples, taken 
from Ref. 4, where possible exceptions have also been pointed 
out. Moreover, we have shown that the second derivative of 
the fidelity is a useful quantity for the quantitative analysis 
of the scaling properties of the system at the transition point. 
From the finite size scaling behavior of S{g) we have indeed 
estimated the critical exponent for the correlation length, find- 
ing agreement with the explicit calculation. Finally, for one of 
the considered examples, we have analyzed both concurrence 
and single site entanglement. Although the latter quantities 
did provide signatures of the undergoing quantum phase tran- 
sitions, this information was hidden in high order derivatives, 
making it much more difficult to extract than the fidelity (or 
its second derivative). 

The fidelity analysis has the advantage of providing a uni- 
fied framework to detect very different types of phase transi- 
tions. For the concrete examples analyzed here, singularities 
in the fidelity are related to level crossings in the transfer op- 
erator, recovering the known transition mechanism for matrix 
product states. The results presented here further contribute to 
demonstrate the generality of the fidelity approach toguantum 
phase transitionsii, supporting the findings of Refs.llB 
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" By using Uhlmann's fidelity for mixed states, one can extend this 
kind of analysis to finite temperature systems and therefore even 
to classical phase transitions. See P. Giorda and P. Zanardi, in 
preparation. 

Considering the fidelity \{g — 5\g + 5) \ corresponds to taking a 
variation 25 between the states. This gives rise to the factor 4 in 
formula ^ and differs from Refs.QQ. 

The correlation length in MPSs is usually given by the formula 
= 1/ In IA1/A2I, where Ai (A2) is the eigenvalue of Et{g) of 
(second) largest modulus. In this case, however, the degeneracy of 
A2 = A3 = — 1 does not allow to apply this expression. 
Note that |0)®^ and [^(ffc))®^ become orthogonal in the TDL, 
in accordance with the normalization factor l/\/2. 
In the limit -> 00 one has {X{g)) ^ for 3 G [0,g'c), 
^ {gTI[l + {g'.f% and {X{g)) 2<7"/(l + g^^) 
for g'^[0, gc]- Hence, for example, limg_^Q- dg{X) = 2n!, while 
lim,_o+ d^{X) = 0. 



